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We analyze in detail recent experiments on ultracold dilute 87 Rb- 40 K mixtures in Hamburg and 
in Florence within a mean-field theory. To this end we determine how the stationary bosonic and 
fermionic density profiles in this mixture depend in the Thomas- Fermi limit on the respective particle 
numbers. Furthermore, we investigate how the observed stability of the Bose-Fermi mixture with 
respect to collapse is crucially related to the value of the interspecies s-wave scattering length. 
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I. INTRODUCTION 



Six years after the first experimental achievement of Bosc-Einstein condensation (BEC) of trapped atomic gases in 
1995 fermionic atomic gases were brought together with bosonic atoms to quantum degeneracy in a 7 Li- 6 Li mixture 
[1, 2], 23 Na- 6 Li mixture [3], and 87 Rb- 40 K mixture [4]. In contrast to a pure Bose system, quantum degeneracy in 
a Fermi system with only one spin component means for T <C Tp = Ep/ks that all energy states below the Fermi 
energy Ep are occupied with one fermion each, whereas all states above Ep remain empty. The main problem to 
achieve quantum degeneracy in a Fermi gas is the inability of fermions to be directly evaporatively cooled. Due to the 
Pauli exclusion principle, fermions in the same spin polarized hyperfine state are not allowed to be close together, so 
that they can not collide via short-range contact interaction to rethcrmalize the gas during the evaporative cooling. 
This handicap was circumvented in the experiment of Ref. [5], where a mixture with two different spin states of 40 K 
was simultaneously evaporated by mutual cooling. In combination with a Bose gas the fermions are sympathetically 
cooled by elastic interactions with the bosons in the overlapping region [4, 6]. 

Beside the exploration of quantum degeneracy, one is also interested in studying how the two-particle interaction 
influences the system properties. Mixtures with a strong interspecies interaction, with the prominent example of 
4 He- 3 He liquid [7, 8], lead to new phenomena like phase separation or BEC-induced interactions between fermions. 
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Depending on the nature of the interspecies interaction, a repulsion between bosons and fermions tends to a demixing 
in order to minimize the overlapping region [9], whereas in the case of an attraction the mixture can collapse as 
long as the particle numbers are sufficiently large [6, 10]. The possibility of superfluidity in a Fermi gas, especially 
the predicted BEC-BCS crossover between BCS-typc superfluidity of Cooper pairs of fermionic atoms and BEC of 
molecules was recently probed [11, 12]. A Feshbach resonance was used to tune the interaction strength between 
fermionic atoms of two different spin states, characterized by the s-wave scattering length a, from effectively repulsive 
(a > 0) to attractive (a < 0). On the BEC side (a > 0) of the magnetic field resonance, there exists a weakly bound 
molecular state. Fermionic atoms are bound into bosonic molecules which can condense at sufficient low temperatures. 
On the BCS side (a < 0) of the resonance two fermions with different spin states form a loosely bounded Cooper 
pair. The condensation of fermionic atom pairs was observed on both the BEC and the BCS side of the Feshbach 
resonance in the experiment [11, 12]. Furthermore, the system properties were observed to vary smoothly in the BEC- 
BCS crossover regime. An alternative and complementary access to Fermi superfluidity is expected from quantum 
degenerate Bose- Fermi mixtures where an effective interaction between fermions is mediated by the bosons [13, 14], 
similarly to phonons in a solid-state superconductor. 

Another recent and fast growing field is the investigation of ultracold boson-fermion mixtures trapped in an optical 
lattice which is created by standing waves of the electric field of counterpropagating laser beams. The atoms can be 
confined to different lattice sites and, by varying the laser intensity, their tunneling to neighboring sites as well as 
the strength of their on-site repulsive interactions can be controlled. In the case of a pure ultracold Bose-Einstein 
condensate with repulsive interaction, held in a three-dimensional optical lattice potential, a quantum phase transition 
from a superfluid to a Mott insulator phase was observed as the depth of the lattice is increased leading to a suppression 
of the tunneling between neighboring lattice sites [15]. The presence of fermionic atoms together with the Bose-Einstein 
condensate makes the system more complex and richer in its behavior at low temperatures. It has been predicted that 
novel quantum phases in the strong-coupling regime occur which involve the pairing of fermions with one or more 
bosons or bosonic holes, respectively, when the boson-fermion interaction is attractive or repulsive [16]. Depending on 
the physical parameters of the system, these composite fermions may appear as a normal Fermi liquid, a density wave, 
a superfluid, or an insulator with fermionic domains. Instead of varying the lattice potential depth the transition 
from a superfluid to a Mott insulator in bosonic 87 Rb can be shifted towards larger lattice depth by adding fermionic 
40 K which interacts attractively with rubidium and therefore increases the effective lattice depth [17, 18]. 

Ultracold trapped boson-fermion mixtures were investigated with respect to a demixing of the components [19- 
21] and to a collapse due to the interspecies attraction [20-24]. Furthermore, the time-dependent dynamics of the 
collapse [25] and finite temperature effects on the stability in a boson-fermion mixture were also studied [26]. Our 
theoretical investigation is based on the experiments with a 87 Rb- 40 K boson-fermion mixture in a harmonic trap, 
which were performed in Hamburg [6] and in Florence [4, 10]. After cooling down at a temperature below 1 ^iK, a 
condensate of 10 6 (2 • 10 5 ) 87 Rb atoms coexisting with 7.5 • 10 5 (3 • 10 4 ) 40 K atoms with a quantum degeneracy of about 



Hamburg Experiment Florence Experiment 



m B = 14.43 • 1(T 26 kg 



m F = 6.636 • 10" 



kg 



(1BF 



a BB = (5.238 ±0.002) nm 
15.0 nm a BF = (-20.9 ±0.8) nm 



mass of S7 Rb atom 
mass of 40 K atom 
s-wave scattering length 
(bosons <-> bosons) 
s-wave scattering length 
(bosons <-> fermions) 

radial trap frequency (bosons) ui B ,r = 2tv ■ 257 Hz uJ B , r = 2tv ■ 215 Hz 

axial trap frequency (bosons) lo b , z = 2ir ■ 11.3 Hz lo b , z = 2-7T • 16.3 Hz 

radial trap frequency (fermions) u)f,t = 2tv ■ 379 Hz lof, t = 2rr • 317 Hz 

axial trap frequency (fermions) lof, z = 2tt ■ 16.7 Hz lof, z = 2tt ■ 24.0 Hz 

number of bosons N B = 10 6 N B = 2 ■ 10 5 



number of fermions 



N F = 7.5 • 10 5 



N F = 3 • 10 4 



TABLE I: List of parameters of the experiments with a 87 Rb- 40 K boson-fermion mixture. The values are taken from 
the experiments in Hamburg [6] and in Florence [4, 10, 27]. 



T/Tp = 0.1 (T/Tp = 0.3) was achieved in the Hamburg (Florence) experiment. The parameters of both experiments 
are summarized in Table I. The distinct values for the interspecies s-wave scattering length as f employed to describe 
each experiment are worth a detailed explanation since this parameter is of great importance for the system, especially 
for the stability of the mixture against collapsing. An overview of different values for clef and their determination 
method along with a reference are shown in Table II. A comparison of the incompatible values for asF shows the 
need of further investigation in this field. 

Now we give an outline of the contents of the paper. At first, we develop in Section II a mean-field theory for 
an ultracold dilute gaseous boson-fermion mixture within the functional integral approach to many-body theory. 
By splitting the Bose fields into background fields and fluctuation fields, we derive an effective action of the Bose 
subsystem without integrating out the fermionic degrees of freedom. Its extremization at zero temperature yields two 
coupled equations of motion, one for the condensate wave function and another one for the Green function of the 
fermions. By evaluating the fermionic Green function for a stationary BEC within the semiclassical approximation, 
we obtain how the time-independent Gross-Pitaevskii equation for the condensate wave function is modified in such 
a mixture. 

In Section III we apply the Thomas-Fermi approximation by neglecting the kinetic energy of the bosons and 
obtain an algebraic Gross-Pitaevskii equation. With the help of its solution we determine the density profiles of both 
components in a 87 Rb- 40 K mixture, where the contact interaction is repulsive between the bosons and attractive 
between both components. 

Finally, we investigate in Section IV the stability of the Bose-Fermi mixture with respect to collapse by numerically 
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a B F /ciBohr Method of determination Reference (year) 

measurement of the elastic cross section for collisions between 

41 K and 87 Rb in different temperature regimes and following [28] (2002) 

mass scaling to the fermionic 40 K isotope 

, ,„ measurement of the rethermalization time in the mixture in 
-330±i^ M ( 2002 ) 

[4, 10] after a selectively heating of 87 Rb 

measurement of the damping of the relative oscillations of 40 K 
— 410±|i [10] (2002) 

and 87 Rb in a magnetic trap 

-395 ± 15 mean-field analysis of the stability of the mixture in [4, 10] [27] (2003) 

magnetic Feshbach spectroscopy of an ultracold mixture of 
-281 ±15 [29] (2004) 

40 K and 87 Rb atoms 

cross dimensional thermal relaxation in a mixture of 40 K and 
250 ± 30 87 Rb atoms after a increase of the radial confinement of the [30] (2004) 
magnetic trap, here only \a,BF/ao\ 

mean- field analysis of the stability, based on [24] , of the mix- 
-284 [6] (2006) 

ture in [6] 

extensive magnetic Feshbach spectroscopy of an ultracold mix- 
-205 ±5 [31] (2006) 

ture of 40 K and 87 Rb atoms 



TABLE II: List of several published values of the s-wave scattering length between 87 Rb and 40 K including their 
determination method and their reference. 

evaluating the effective action for a trial Gaussian density profile of the condensate. We compare our results, which 
strongly depend on the value of the Bose-Fermi s-wave scattering length, with the experiments on 87 Rb- 40 K mixtures 
in Hamburg and Florence. The good agreement of our theoretical results with the measurement in the Florence 
experiment allows us to fit <jbf which describes the Hamburg data quite well, but remains incompatible with the 
value used in Florence. 



II. DERIVATION OF GROSS-PITAEVSKII EQUATION 



In this section we summarize briefly the mean-field theory of boson-fermion mixtures [14, 24] and derive within 
the functional integral formalism a coupled set of differential equations for the condensate wave function and the 
fermionic Green function. 
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A. Grand-Canonical Partition Function 

We consider a dilute gaseous mixture of ultracold bosonic and fermionic atoms. In order to obtain quantum 
statistical quantities for such a Bose-Fermi mixture, we use the grand-canonical partition function in the functional 
integral formalism. Thus, we integrate over all possible Bose fields Vb( x , r), ip B (x,r) and Fermi fields ip F (x,T), 
V>f(x, t), which are weighted by a Boltzmann factor [14, 32, 33]: 



z = <j> Vip* B <j> v^ B <j> Vip* F <j> vi> F e -A[r B ^B,r F ^F]/h_ 



(1) 



The complex fields %(j b (x,t), tp B (x, r) represent the bosons and are periodic on the imaginary time interval [0, h/3], 
whereas the fermions are described by Grassmann fields ip F (x,T), V"f(x, r) which are antipcriodic on this interval: 

V>*(x,?i/?)=eV*(x,0), V(x,fi/3) =eV(x,0). (2) 

Here e = ±1 holds for bosons and fermions, respectively. The total euclidean action of a Bose-Fermi mixture consists 
of three parts: 

A^b^b^f^f] = Ab[^b,^b}+ A f [iI> f ,^f]+ AbfWb^b^f^f]- (3) 
The first term describes the bosonic component of the mixture: 



rnp r- 

AbWb^b] = j dr J d 3 x^* B (x,T) 



^b(x,t). (4) 



It contains the Legendre transform, the kinetic energy, the external trap potential Vb(x), the chemical potential fx B 
to fix the boson number, and the strength g B B — 4:Trh 2 a B B/mB of the contact interaction between two bosons with 
the s-wave scattering length a BB - As we deal with a dilute gas, an interaction between more than two particles 
is negligible, and we can restrict ourselves to the short-range contact interaction. Since the Pauli principle forbids 
fermions in the same hyperfine state to be close together and therefore to collide via contact interaction, we can write 
the corresponding action term for the fermionic component of the mixture as 



rh/3 ,■ 

AfW f Af\= J dr J d 3 xxP* F (x,T) 



2 



^-2^ A + ^ (X) -^ 



^f(x,t), (5) 



where Vf(x) and [If denote the external trap potential and the chemical potential for fermions. As we assume a 
situation, where the bosonic and fermionic atoms cannot be transformed into each other, each of both species has its 
own chemical potential. The last term of the euclidean action (3) 

/■ft/3 r 

AbfI^b^bAfAf] = gBF J dr J d 3 x\tp B (yi, r)| 2 |V>f(x, t)| 2 (6) 

describes the contact interaction between bosons and fermions, where its strength gBF is related to the s-wave 
scattering lengths a B p via [34, 35] 

2 m B+ m F , > 

g B F = 2nn,a B F • (7) 

tub rriF 
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B. Background Method 

In order to account for the fact that the bosons in the mixture can condense, we apply the background method of 
field theory [33, 36, 37] and split the bosonic fields tjj B (x,T), ip B (x,r) into two parts: 

i/>% (x, r) = f B (x, r) + 6ip* B (x, r) , Vb (x, t) = * B (x, r) + 6ip B (x, r). (8) 

The first part represents the background fields *f? B (x., r), ^(x, r) whose absolute square is identified with the density 
of the condensed bosons. The second part consists of the fluctuation fields 5tp B (x, r) , #V-b( x > t ) of the Bose gas 
describing the excited bosons, which arc not in the ground state. Using the decomposition (8), we expand the 
cuclidean action (3) in a functional Taylor series with respect to the Bose fields i(> b (x,t), Vb( x ; r ) around the 
background fields ^(x, t), "^(x, r). In this work we restrict ourselves to the Gross-Pitaevskii theory, i.e. we consider 
the euclidean action A only up to the zeroth order in the fluctuation fields 5ip* B (x, r) , Sips (x, r) which is equivalent 
to evaluate the euclidean action (3) at the background fields: ■4[\I , £ + Sip B , + StpB, ipp, iPf] ~ -^[^b; ^b, i^p, iPf]- 
Thus, the bosonic functional integration in Eq. (1), whose integration measure transforms according Dip B (x, r) = 
VSip B (x, t) and £>V-b( x , t) = VStpsi^-, t), can be dropped. In this way we obtain for the effective action r[*g, ^ B ] = 
-\nZ[^* B ,^ B }/P the result 

r[^, *b] = ^ *b] - ^ lnZ F [^, * B ], (9) 

where 

Z F [** B ,* B ] = <j,Vip* F j > T>^ F e- AF ' et ^' <5!B ' r ^ F ^ h (10) 

represents the functional integral over the Fermi fields resulting in a pure functional of the Bose background fields 
^ b (x,t), fg(x,r). The effective cuclidean action Af^^s-i^Bi^p^f] depending on the Fermi fields i/)^(x,r), 
■0f(x, t) is summarized by 



,ft/3 , 

^[^,^,#,^1 = y dry d 3 x^(x,T) 
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fi-TT + i?F,cff (x, t) - /i F 



Vf(x,t), (11) 



where the effective one-particle Hamilton operator for fermions reads: 

h 2 

Hp, eS ( X ,r) = -- A + V F (x)+. 9BF |* B (x,r)| 2 . (12) 

imp 

C. Coupled Equations of Motion 

The grand-canonical free energy T is the extremum of the effective action r[5'|j, ^ b\ with respect to the background 
fields **(x, r), W(x, t). Using Eqs. (4) and (9)-(12), this yields the time-dependent Gross-Pitaevskii equation: 

d h 2 



h^- - — — — A + Vb(x) - /i B + |*b(x, t)| 2 - g BF G F (x,T;x,T) 
or 2m,B 



*b(x,t) = 0, (13) 
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where the fermionic two-point function is defined as follows: 

Gf(x, t; x', t') = — j> Vr F j> V^ F ^(x, r) ^(x', r') e -^.^[*i.*B,^,^]/ft (14) 

This imaginary time-dependent nonlinear Schrodinger equation (13) of the condensate wave function ^(x, r) repre- 
sents a partial differential equation, where the nonlinear terms are due to both interactions. The last term results from 
the boson-fermion interaction, whereas the other terms have the conventional Gross-Pitaevskii form for a condensate 
[38, 39]. The fermionic two-point function is at the same time a Green function for the fermion fields, thus, it obeys 
the linear inhomogeneous Schrodinger equation for fermions 

G f (x,t;x',t') = TiS(x-x')S^(t-t'), (15) 

where the inhomogeneity consists of the delta function in space and the antiperiodic repetitive delta function in 
imaginary time 

oo 

5^{t-t')= {-l) n 5{T -t' + nh(3). (16) 

n— — oo 

The Gross-Pitaevskii equation (13) for the condensate and the inhomogeneous Schrodinger equation (15) for fermions 
form a set of coupled equations in imaginary time. Therein the condensate wave function ^s(x, r) in the Gross- 
Pitaevskii equation (13) is modified by the fermionic Green function G_f(x, t; x, r) and, vice versa, the condensate 
wave function \I/b(x, t) influences the fermionic Green function Gf(x,r;x',r') in the inhomogeneous Schrodinger 
equation (15). A Wick rotation r = it and the omission of the chemical potentials in Eqs. (13) and (15) leads to two 
coupled equations of motion in real time, which describe the dynamics in the Bose-Fermi mixture. 



dr 



#F,eff(x, T) - [ip 



D. Fermionic Green Function 



In the following we obtain an approximative solution for the fermionic Green function Gj?(x, t;x',t')- To this 
end we restrict ourselves to a stationary BEC so that the imaginary time-dependent Gross-Pitaevskii equation (13) 
reduces to 

h 2 



2m,B 



A + Vb(x) - fi B + gss|*s(x)| - gs F G p(pc, t; x, t) 



*b(x) = 0, 



and the effective fermionic Hamiltonian (12) no longer depends explicitly on the imaginary time r: 

F.cff (x) = - A- A + V F (x) + g BF | * B (x) | 2 . 

zm F 



(17) 



(18) 



Now we assume that both the fermionic trap potential Vf(x) and the condensate wave function ^E'b(x) have only a 
weak spatial dependence so that the semiclassical approximation can be applied. Therein the eigenvalue problem 



Hp >eff (x) Vk(x) = £(k, x) Vk(x) 



(19) 



is approximately solved by plane waves 



^k(x) = 



, ? k x 



(2tt) 3 / 2 



and the scmiclassical energy spectrum 



ft 2 k 2 



£7(k,x) = ^ + V F (x) + 3 bf|*b(x)| 2 . 

Thus, the semiclassical solution of Eq. (15) yields for the fermionic Green function 
d 3 k 



(20) 



(21) 



G ,(x,.;xV)=/^- 



»k.(x-x') 



0( - T _ r /) e -[B(k,(x+x')/2)-MP](r-r'-ft0/2)/ft _ Q(y _ r ) e -[J5(k,(x+x')/2)-^](r-r'+ft/3/2)/ft 

X 2cosh/3[i;(k, (x + x')/2) - /z f ]/2 ' ^ 

The limit of equal imaginary times follows from r' j r: 



Gf(x,t;x',t) = - 



= ik-(x-x') 



(2?r) 3 e /3[S(k,(x+x')/2)- A i J ,] + ! 

Thus, the fermionic particle density fij?(x) = -Gp(x, r; x, r) is given by 

d 3 k 1 



n F (x) 



(2tt) 3 e /?[fi 2 k 2 /2m F -AF(x)] + ! 



with the local chemical potential 



Mf(x) = Mf - V>(x) - |*b(x)| 5 



(23) 



(24) 



(25) 



which represents the kinetic energy of the fermion in the highest energetic state, when it is located at the space point x. 
In order to evaluate the integral (24) in momentum space, we apply spherical coordinates and the substitution 
e(k) = H 2 k 2 /2m F : 

1 / m F \ 3 / 2 [°° dee 1 ' 2 



n F (x) = 



r(3/2) \2nh 2 J J e^-^Wl + l' 
In the low-temperature limit T J, 0, the Sommerfeld expansion [40] of Eq. (26) yields in the lowest order 

n F (x) = (t9(/i F (x)) Af 2 (x) 



(26) 



(27) 



with the abbreviation k = (2m F ) 3 ' 2 /6-K 2 h 3 . Finally, we insert the result for the fermionic particle density (27) into 
the Gross-Pitaevskii equation (17): 

h 2 



-A + Vb(x) - MB + .g Bi3 |1< B (x)| 2 + ng BF 6 (/*f( x )) MfV) 



*b(x) = 0. 



2tob 

Note that this stationary Gross-Pitaevskii equation (28) follows from extremizing the effective action 

h 2 



(28) 



r[* s ,* s ] = | rf 3 x|^( x ) 



■^r-A + y s( x)-M S + ^l^(x)| 2 



* b (x)-^9(Af(x))4 /2 (x)| (29) 



with respect to ^* B (x). This effective action can also be obtained from Eqs. (4) and (9)-(12) by evaluating the 
fermionic functional integral (10) in the semiclassical approximation for a stationary BEC and by performing the 
Sommcrfcld expansion in the low-temperature limit. Identifying the extremum of the effective action T with the 
grand-canonical free energy T, the number of bosons and fermions are obtained from (29) via the normalization 
condition 



N j = -^=JdPxn j ( ? 



j = B,F, (30) 



where the particle densities of bosons and fermions read 



n B (x) - |*s(x)| 2 , (31) 
n F (x) = k9(/2 f (x)) 4 /2 (x). (32) 



III. DENSITY PROFILES 



In this section we solve the stationary Gross-Pitaevskii equation (28) in the Thomas-Fermi approximation in order 
to calculate the boson and the fermion density distribution for the parameters of the experiment in Hamburg [6]. 
There the 87 Rb- 40 K boson-fermion mixture is confined in a three-dimensional rotationally symmetric harmonic trap 

^(x) = ^i(^ r r 2 +a;^ 2 ), i = B,F, (33) 

where we use cylindrical coordinates {r, (j), z}. The frequencies for both species are related by UF,k — \frnB~fmF uJs.k 
for k — r,z, so that tub uj% k /2 = mp top k /2. 



A. Thomas-Fermi approximation 

Assuming that the potential and interaction energy are larger than the kinetic energy, we can use the Thomas-Fermi 
approximation, where the kinetic term in the Gross-Pitaevskii equation can be neglected. With this approximation 
the Gross-Pitaevskii equation (28) reduces with the help of Eqs. (31) and (32) to an algebraic equation with respect 
to the bosonic particle density ris(x) [41]: 

Vb (x) - [is + 9bb n B (x) + g B F n F (x) = 0. (34) 

The last term in Eq. (34) stands for the fermionic particle density 

n F (x) = k8 (fi F - Vf(x) - 9bf n B ( x )) [fJ-F - Vp(x) - g B F n B (x)] 3/2 , (35) 

which modulates the bosonic density profile ns(x) and vice versa. 
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B. Vanishing Boson-Fermion Interaction 



Now we discuss the special case of vanishing boson-fermion interaction gsF — * 0, where we mark all quantities with 
the upper index (0) . On the one hand, Eq. (34) reduces then to the well-known particle density of a pure BEC in the 
Thomas-Fermi approximation 



»£>(*) 



1 e (/4 0) - vb(x) 



9bb 



(0) 



Vb(x) 



(36) 



provided that gas > 0. In case of an attractive interacting BEC with negative gss-, the third term in Eq. (34) cannot 
be balanced in the trap center at x = by the remaining chemical potential ji^ leading to a collapse of the BEC. 
The particle density (35) of fermions, on the other hand, becomes independent of that of bosons: 



,(0) 



(x) = «e(74 0) 



Vf(x) 



v>(x) 



3/2 



(37) 



Thus, the clouds of the BEC and the Fermi gas coexist undisturbed from each other. Setting the particle densities 
(36) and (37) zero yields the Thomas-Fermi-radii of the BEC and the Fermi cloud along the respective axes: 



d(0) _ 



\ 



2fi 



(0) 



i = B,F. 



(38) 



The chemical potentials (J^ and are determined by the particle numbers Nb and Np via the normalization (30): 



(o) 
Mb 



15cibbN b \ 2/5 hu} B 
L R 



(39) 



where Lj = ^Jh/rrn u>i denotes the geometrical average of the oscillator lengths L^k — y/h/rrii w^k. Inserting the 
chemical potentials (39) into the respective particle densities (36) and (37), the maximum of the latter in the trap 
center results in 



n<°>(0) 



(15iVB) 2/5 

„ 3/5 ;12/5 ' 

8to bb l b 



nf{0) 



2N 



1/2 



(40) 



In the Thomas-Fermi approximation the BEC needs a non- vanishing repulsion between the bosons in order to prevent 
the density from becoming infinite or complex, which indicates a collapse of the BEC, whereas the quantum pressure 
in the noninteracting Fermi gas preserves the latter from a collapse. The maximum fermionic particle density increases 
faster with the particle number than the maximum bosonic particle density, whereas the latter depends not only on 
the oscillator length but also on clbb> the parameter for the bosonic interaction strength. The particle densities of 
both species in an undisturbed BEC and Fermi gas are plotted for typical particle numbers Nb and Np as thin lines 
in Figure 1. 



C. Non- Vanishing Boson-Fermion Interaction 

After having treated the case gsF = in the previous section, we now discuss the algebraic equation (34) for a 
BEC which interacts with fermions. In order to plot the corresponding bosonic and fermionic densities in Figure 1, 
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1.5 - 



0.5 - 




9bb n B (x)/^ 
\g B F\ n F {yi)lix% ) 

(0)/ \ / (0) 

9bb ij, W//iJ, j 



\9bf\ »f ( x )/Vb 
N B = 1.74 • 10 6 
N F = 7.83 • 10 6 



1.5 



FIG. 1: Comparison of the dimensionless particle densities for bosons and for fermions between a disturbed mixture (thick 
lines) and an undisturbed BEC and Fermi gas (thin lines). The densities are plotted versus the coordinates r at the plane 
2 = and 2 at the line r = 0, respectively, in units of the Thomas-Fermi radii of the undisturbed BEC. The densities of the 
disturbed mixture are the solution of Eqs. (34) and (35), whereas those of an undisturbed BEC and Fermi gas are given by 
Eqs. (36) and (37), respectively, for an example of typical particle numbers N B , Nf of the Hamburg experiment [6]. 



we insert the fermionic density (35) into the algebraic equation (34) and obtain: 

Vb(x) - fi B + 9bb n B (x) = -g B F k \p F - V F (x) - g B F ns(x)] 3/2 . 
Squaring (41) leads to a cubic equation with respect to the condensate density n_e(x): 

n|(x) 



(41) 



g% B 3\mf-V f (x)]\ 2 f 2g BB [V B (xL) - » B ] 3[ Mf - V F (x)] 2 
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It is of the form 



ax 



bx + c = 



(43) 



with real coefficients a, b, c and its three solutions x\, X2, £3 can be found with the help of the Cardanian formula [42]. 
The first solution x\ remains always real, whereas the other solutions xi and X3 are either real for 



a 
~9 



n 3 



< 



a 
27 



ab c 
¥ ~ 2 



(44) 



or conjugate complex otherwise. Although all three solutions obey the cubic equation (43), i.e. Eq. (42), only the 
last solution x% turns out to satisfy the algebraic equation (41). Inserting the two other solutions into the algebraic 
equation yields the right value but opposite signs on both sides of the algebraic equation (41). The reason is that 
Eq. (41) represents the root of the cubic equation (42), therefore an equation of the order 3/2, which reduces the 
number of possible solutions. 

The resulting particle densities for both species are plotted as thick lines in Figure 1 for typical particle numbers N B 
and Np using the experimental parameters in Table I. In order to determine the solution outside the BEC, especially 
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the particle density of fcrmions, we must not forget the origin of the algebraic equation (34) from the Gross-Pitaevskii 
equation (28). Although the solution n_e(x) becomes negative outside the BEC, we have to set it equal zero there, as 
the density must be positive by definition. Thus, the condensate wave function \&b(x) becomes also zero outside the 
BEC and fulfills the Gross-Pitaevskii equation (28) in a trivial way regardless of the expression inside the brackets, 
which is, apart from the neglected kinetic energy term, the algebraic equation (34). Hence, the fermionic particle 
density outside the BEC is described by Eq. (35) with ns(x) =0 without obeying the algebraic equation (34). 
The particle density of the BEC in the Thomas-Fermi approximation can be written in the following form: 

n B (x) = — !— [hb - Vb(x) - g BF n F (x)] . (45) 

9BB 

From the fermionic particle density (35) we deduce that its maximum occurs in the trap center x = due to the 
negative gBF- Thus, the BEC possesses its largest density also in the trap center. In other words, the particle 
densities of both species intensify each other in the overlapping region due to the strong attraction between bosons 
and fermions as is shown in Figure 1. With increasing distance from the trap center both densities (35) and (45) 
decrease quickly within their overlap due to the interspecies interaction terms and due to the increasing trap potentials, 
whereas outside of the overlap only the latter reason is responsible for their decreasing. This behavior is shown in 
Figure 1 for the Hamburg experiment. Which of both clouds has the larger extension depends on the particle numbers 
Nb and N F and therefore on the chemical potentials \xb and [i F . Usually, the BEC is surrounded by the Fermi gas 
unless Nb 3> N f . Both chemical potentials, which represent the total energy of a particle of the corresponding 
species, are smaller compared with /j,^ and fjp^ of the undisturbed BEC and Fermi gas, since the particles possess, 
besides the neglected kinetic energy, the potential energy due to the trap and the intraspecies interaction energy an 
additional negative interaction energy due to the interspecies interaction. Figure 2 shows these energies of a boson 
in units of V B (RB,r,0) = V B (0, Rb.z) = I^b f° r the Hamburg experiment, where we see that hb ~ 0.6 These 
reduced chemical potentials lead to decreasing Thomas-Fermi radii for both species. Indeed, one can see in Figure 
1 a reduction of the Thomas-Fermi radii of the BEC to Rs.k ~ 0-SRg\. Hence the attractive interaction between 
both species leads to an additional confinement of both the BEC and the Fermi gas within their overlap. 

A direct quantitative comparison with the density of the experimental probe is not possible since a measurement of 
the density in the trap fails due to the smallness of the probe. Only after the expansion of the probe, when the trap 
potential was suddenly switched off, its density can be measured taking the absorption image of the optical density of 
87 Rb and 40 K. This would allow to compare the densities only qualitatively, since the density in the probe is changed 
due to both the ballistic expansion and the dynamics in the probe. 



13 




FIG. 2: Comparison of the kinetic energy Ekin = -S 2 A$B(x)/(2ms *1>b(x)) of a boson, with its remaining energies, namely the 
intraspecies interaction energy £int,BB = <7bb ns(x), the interspecies interaction energy -Eint.BF = |<7bf| tif(x), the potential 
energy E pot = Vs(x) due to the trap and the chemical potential -Echom = fJ-B- All energies are related to the value of the 
potential energy at the boundary of the undisturbed BEC cloud and are plotted versus the coordinates r and z within the BEC 
cloud for the situation in Figure 1. 



D. Validity of Thomas- Fermi approximation 



In order to check the validity of the Thomas-Fermi approximation, we have plotted in Figure 2 the kinetic energy 
of a boson 

E ^ ~ 2m s v M x) (46) 

in units of the value ^ of the potential energy at the boundary of the undisturbed BEC. The kinetic energy is, 
indeed, negligible in a wide bulk range from the trap center to just before the boundary of the disturbed BEC. Thus, 
the Thomas-Fermi approximation gives very accurate results except in the outermost 10% of the Thomas-Fermi radii. 
The kinetic energy diverges at the BEC boundary because the condensate wave function ^(x), which represents 
the square root of the particle density ns (x) , occurs due to the Laplacian derivative in Eq. (46) in the denominator, 
which becomes zero at the BEC boundary. It is obvious that the proper solution of the Gross-Pitaevskii equation 
(28) would match the one of the algebraic equation (34) from the trap center to just before the BEC boundary, where 
it tends smoothly to zero, thus improving the sharp bend in the graph of the fermionic particle density there to be 
smooth. 



E. Complex solutions 



As already discussed in Subsection IIIC, the solution n_e(x) of the algebraic equation (34) does not always remain 
real by varying the chemical potentials fis and fip. From the condition (44) we extract a straight line in the 
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FIG. 3: (a) Stability diagram with respect to the chemical potentials. The solid line separates the stable region on the left from 
the unstable region on the right according to Eq. (47) . The pictures (b)-(d) show a sequence how the bosonic and fermionic 
particle densities versus the coordinates r at the plane z — and z at the line r = 0, respectively, change on the road from 
stability to instability. The {^b, n F )-pairs belonging to (b)-(d) are located on the dashed line in (a) and are equally spaced 
A/ub = 5 • 1CP 32 J apart giving rise to particle number differences of about ANb ~ 1.8 • 10 5 and AJVf w 3 • 10 4 , respectively. 
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which separates the half plane with a complex solution from that with a real solution. This situation is depicted 
for the Hamburg experiment in Figure 3 (a). Due to the results of the subsequent Section IV, we can associate the 
complex solution to an instability of the mixture with respect to collapse. Figures 3 (b)-(d) demonstrate the change 
of the density profiles ns(x) and u,f(x) on the road from stability to instability by increasing the particle numbers 
with ANb ~ 1-8 ■ 10 5 and ANp ~ 3 • 10 4 each. The respective pairs (hb, Hf) lie on the dashed line, which is arranged 
perpendicular to the critical solid line in Figure 3 (a), so for Figure 3 (b) just below the critical line, for Figure 3 
(c) on, and for Figure 3 (d) just above. Figure 3 (b) shows a stable configuration, Figure 3 (c) is on the boundary 
between stable and unstable where the densities pile up at the trap center to a peak, and in Figure 3 (d) the real 
part of the densities is chopped off at the trap center, so that there an imaginary part of both densities occurs. This 
imaginary part starts to appear simultaneously for both components at the density maximum x = and grows in 
magnitude and in extension out of the trap center with increasing Nb and Np. Such an anomaly in the densities was 
observed in the Hamburg experiment, where the evolution of an overcritical mixture is shown in Figure 1 of Ref. [6]. 



Here we determine the stability border both within the Thomas-Fermi approximation and, in a separate variational 
calculation, beyond the Thomas-Fermi approximation. The stability border turns out to depend strongly on the 
value of the interspecies s-wave scattering length. Therefore, comparing our theoretical results with the experimental 
measurements allows to extract a trustworthy value for this crucial s-wave scattering length. 



In Subsection HIE we described in detail the behavior of the particle densities ns(x) and rip(x) as the solution of 
the stationary Gross-Pitaevskii equation within the Thomas-Fermi approximation for varying particle numbers Nb 
and Np . We assigned the emergence of a complex density to a loss of the stability against collapse and found a border 
in form of a line in the (^b, MF)-planc, which separates the stable and unstable regions as shown in Figure 3 (a). In 
order to obtain a stability diagram in the (Nb, Np) -plane, we have evaluated the corresponding particle numbers from 
the chemical potentials by integrating out the respective particle densities according to the normalization condition 
(30). The result for the Hamburg and the Florence experiment is given by the dot-dashed line in Figures 4 and 
5, respectively. Mixtures with particle number pairs (Nb,Nf) below this line are stable whereas particle number 
pairs above the line indicate an unstable mixture tending to collapse. The critical particle numbers of both species 
behave, roughly spoken, inversely proportional to each other in a wide range. Whereas the critical number of bosons 
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A. Thomas-Fermi Approximation 
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FIG. 4: Stability diagram for the 87 Rb- 40 K mixture of the Hamburg experiment. The particle number pairs (Nb, Nf) below a 
certain line belong to a stable mixture wheras those above the line indicate an unstable mixture tending to collapse. The thin 
lines correspond to the quantum mechanical limit with the ratio Aqm = u) z /lu t and the thick lines represent the Thomas- Fermi 
limit with the ratio Atp = {uj z /uj r ) 2 ■ The points are obtained in the experiment by analyzing decay series in various particle 
number regimes and are assigned to stable mixtures (dots) and to unstable mixtures (crosses) [6] . 



NbciH tends to zero when the number of fermions increases, the critical number of fermions Npait remains finite and 
constant when the number of bosons is enlarged. This situation, where the line in the stability diagram becomes 
vertical, happens when the BEC cloud becomes so large that it surrounds the Fermi gas. This defines a minimal 
number of fermions -/V^min below which the mixture remains stable irrespective of the number of bosons Nb- In 
order to estimate Np :in i n , we have to integrate the fermionic particle density (35) by taking into account the stability 
condition (47) when solving the algebraic equation (41) for the condensate density tib(x). In this way we find that 
the fermionic particle density (35) does not explicitly depend on both chemical potentials \±f and fiB- A numerical 
evaluation for the Hamburg (Florence) experiment yields ATp min = 9,99 • 10 5 (A?ir m i n = 1,08 • 10 4 ) in accordance 
with the dot-dashed instability line shown in Figure 4 (Figure 5). 



B. Variational Method 



Another approach to determine the stability border for a 87 Rb- 40 K mixture is based on extremizing the grand- 
canonical free energy (29) with the local chemical potential (25). In contrast to the Thomas- Fermi approximation in 
Subsection IV A, this approach takes into account the kinetic energy of the condensed atoms. But instead of varying 
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FIG. 5: Stability diagram for the 87 Rb- 40 K mixture of the Florence experiment. The particle number pairs (Nb,Nf) below a 
certain line belong to a stable mixture whereas those above the line indicate an unstable mixture tending to collapse. The thin 
lines correspond to the quantum mechanical limit with the ratio Aqm = u) z /lu t and the thick lines represent the Thomas- Fermi 
limit with the ratio Atp = (u> z /io r ) 2 - The crosses represent mixtures in the experiment which are found very close to the 
instability, whereas the dots indicate stable mixtures [27]. 



the condensate wave function x &s(x) in the grand-canonical free energy (29), which leads to the Gross-Pitaevskii 
equation as described in Section II D, we use the ground-state wave function of a three-dimensional anisotropic 
harmonic oscillator 



*b(x) = 



N B 



^l 2 o?L% 



cxp • 



(48) 



as a test function with variational widths a L B ,k- The dimcnsionlcss factor a, which scales the oscillator lengths of the 
Gaussian function (48), serves as a variational parameter, which extremizes the free energy (29). The test function is 
normalized to N B bosons and obeys for a = 1 the Gross-Pitaevskii equation for a trapped noninteracting BEC: 

h 2 



2m B 



A + Vb(x) - fl B 



*b(x) =0. 



(49) 



Here we assume that the condensate wave function ^s(x) has, also in case of intraspecies and interspecies two- 
particle interactions, qualitatively the shape of a Gaussian curve, as one can read off from Figure 1. Using cylindrical 
coordinates {r, cf), z}, the test function (48) reads 

r 2 + Xz 2 ) 
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Beside a uniform variation of the widths aLs,k by the factor a, we have to consider that the ratio A = (L b , z / 'L b , r ) 2 
could also change due to the interactions. In order to include this, we perform the calculation with two different ratios. 
On the one hand we use Aqm = u z /uj r which stands for the limit of vanishing interactions gsB —* and gsF — > 
and reflects the proper ratio of the oscillator lengths in the ground-state wave function of the quantum-mechanical 
harmonic oscillator (49). On the other hand, we set Atf = (^M) 2 , which represents the Thomas-Fermi limit of 
negligible kinetic energy due to strong intraspecies and interspecies interaction. Inserting the test function (50) into 
the grand canonical free energy (29) reduces the latter from the functional J-[^* B , ^ b] to a function T(a) of the 
parameter a. As the test function (50) is normalized independent of /is, the latter plays no longer a role in T(a) 
and shifts the free energy only by a constant value. On the other hand, the fcrmionic chemical potential fip is needed 
for evaluating the respective fermion number Nf by integrating out the fermionic particle density (35) according to 
the normalization (30). The dependence of the grand-canonical free energy J 7 (a) on the variational parameter a for 
given chemical potentials fis and is shown in Figure 6 for several boson numbers Nb- For Nb < NbctH the free 
energy T(a) possesses a local minimum which corresponds to a metastable state of the mixture. The condensate wave 
function v I / _b(x) has finite equilibrium widths a cq L_B,fc, where a eq denotes the parameter at the local minimum. When 
the boson number exceeds the critical value, i.e. Nb > NbctH, the local minimum disappears so that the widths tend 
to zero in order to minimize T(a). Just this happens when the mixture collapses. Thus, the border between stability 
and instability is given by the condition 

d 2 T(a) 



dT(oe) 

N B = AWit ^ 



da 



da 2 



0, (51) 



i.e. Fict) has a point of inflexion at a — a C rit- The appearance of the local minimum arises from the competition 
between the positive first three terms of the grand-canonical free energy (29) and the negative last term describing 
the influence of the fermions. We have determined the stability border within the variational method in two different 
ways, which we discuss in the following subsections. 



C. Numerical Integration 

The most reliable approach is based on performing the integration in Eq. (29) numerically. We have evaluated 
the critical boson numbers Nscrit for different values of fip in an iterative way until the condition (51) is achieved 
with sufficient accuracy. The result is given by solid lines in the stability diagrams for the Hamburg and Florence 
experiment as shown in Figures 4 and 5, where the resulting particle number pairs (Nb, Nf) are smoothly connected 
with each other. The thick solid line of the Thomas-Fermi ratio Atf = (^z/^r) 2 lies above, but near the dot-dashed 
line from the Gross-Pitaevskii equation in the Thomas- Fermi approximation in Section IV A. This is expected as both 
the thick solid line and the dot-dashed line are evaluated in the Thomas-Fermi limit with the same ratio of the radial 
and axial extension of the BEC cloud. Furthermore, they show the same behavior for very large boson numbers as 
both lines become vertical so that the fermion number remains constant. In the Florence experiment both lines lie 
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FIG. 6: Grand-canonical free energy !F(a)/NBhwB,r versus 
the variational parameter a for various boson numbers Nb- 
A local minimum indicates a metastable state of the mixture. 
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FIG. 7: Comparison of two stability borders whose ratio A 
is 40% smaller and larger, respectively, than Atf in order to 
estimate the proper ratio A which allows the largest possi- 
ble particle number for a stable condensate of the Hamburg 
experiment. 



very close to the crosses of the experiment where the thick solid line fits them better. This is not surprising as the 
s-wave scattering length <jbf in Tabic I is determined with a mean-field analysis in the Florence experiment [27]. The 
thin solid line lies far below the thick solid line and reflects a mixture in the quantum-mechanical limit with the ratio 
Aqm = u z /oj r . This line, which is a good approximation for mixtures with vanishing intraspecies and interspecies 
interactions, is less suitable for the 87 Rb- 40 K mixture as the interaction energy is dominant according to Figure 2. 
Hence, the mixture can be well described in the Thomas-Fermi limit. As the thin solid line does not consider the 
proper ratio A of the oscillator lengths in the test function (50) and, thus, minimizes the free energy less optimally for 
a strong interacting mixture, this line allows much smaller particle numbers in a stable mixture than the thick solid 
line. In order to estimate the ratio A, which yields the stability border with the largest possible numbers for bosons 
and fermions, we have evaluated two lines in the stability diagram for ratios in the neighborhood of Atf, where one 
A is 40% smaller and the other one is 40% larger. These lines for A = 0.6 Atf and A = 1.4 Atf are shown in Figure 7. 
Both lines lie throughout below, but very close to the original one with A = Atf- This indicates that the stability 
border in the Thomas-Fermi limit is stationary at A = Atf where it possesses a maximum. Hence, the ratio Atf turns 
out to be the proper ratio for the 87 Rb- 40 K mixture allowing the largest possible numbers for bosons and fermions 
in a stable mixture. 

D. Expansion with respect to Interspecies Interaction 

Another approach was suggested in Ref. [24], where the local chemical potential (25) in the grand-canonical free 
energy (29) is expanded up to the third order in ggp in order to get rid off the fractional power in the last term of 
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the free energy. This expansion leads to Gaussian integrals with respect to the test function (50): 
**] = / d 3 x ^- |V* B (x)| 2 + y eff (x) |* B (x)| 2 + M |* b(x) |4 + |^ fl(x) |« 



where the terms with respect to the power of |^b(x)| 2 are summarized in the factors 
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(52) 

(53) 
(54) 



Within this approximation it is assumed that the radius of the condensate is much less than the radius of the Fermi 
gas cloud so that the remaining expressions (fj,p — Vf( x ))™/ 2 with n = 1,3,5, arising from the expansion, can be 
expanded in powers of Vf(x)/hp as well. We consider therein only terms which depend on the parameter a. The last 
term in the free energy (52) corresponds to the elastic three-particle collision induced by the interspecies interaction. 
This term with gsF < is responsible for increasing both the bosonic and the fermionic particle density in the trap 
center in order to minimize the free energy. If the central condensate density |^s(0)| 2 becomes large enough due to 
large particle numbers, the positive first three terms in the free energy (52) cannot balance the negative last term 
in order to stabilize the mixture and to prevent it from collapsing. Performing the Gaussian integration in Eq. (52) 
leads to an algebraic equation with respect to the unknown quantities a and Nb- 
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with the factors 
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The condition (51) for the stability border provides two equations allowing to determine both unknown quantities 
for different values of Hf- The result is shown in Figures 4 and 5 by the dashed line for both limits. The thick 
dashed as well as the thin dashed line converges with increasing fermion number Np to the corresponding solid line 
of the numerical integration. But for low Np or, equivalently large Nb, these dashed lines stay below the solid lines, 
where the discrepancy increases with decreasing Np. Moreover, for very small Np the dashed lines show the opposite 
behavior of the solid lines as they tend to zero, which seems to be unphysical. When the radius of the BEC clouds 
increases with decreasing Np, the above mentioned expansion in powers of Vp(x)//j,p, which is done up to the zeroth 
and first order, fails as Vp(x) and \ip become comparable at the BEC cloud boundary. Thus, more orders of the 
expansion are needed to obtain more accurate results. Another reason is that the stability border in the stability 
diagram depends strongly on the interspecies s-wave scattering length cibf according to the scaling law for the critical 
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FIG. 8: Stability border for a new value bbf = ( — 16.82 ± 0.30) nm of the 87 Rb- 40 K mixture in the Hamburg experiment 
[6] . All three lines are obtained by numerical integrations of the fermionic density (35) according to Subsection IV C in the 
Thomas- Fermi limit with the ratio Atf = (cu z /uj r ) 2 . Those stable (unstable) points, which are located above (below) the solid 
line, are equipped with error bars indicating a relative uncertainty of 20% and 30% for Nb and Nf, respectively, according to 
Figure 3 of Ref. [6]. 



numbers of condensate atoms for a fixed ratio between Nb and Np [27, 41]: 

-Wficrit 



(57) 



Thus, such a strong sensitivity of the critical boson number Nscrit with respect to clef and also to gsF due to Eq. (7) 
does not justify an expansion of the local chemical potential /2f(x) with respect to the smallness parameter gsF- 



E. Adjustment of cibf 

The above described strong dependence of iVscrit on clbf allows us to extract a value for asF with great accuracy 
within a mean-field analysis of the stability as performed in this section. Because of the scaling law (57) the relative 
uncertainty of clef amounts to only a twelfth of the relative uncertainty of the critical boson number Nsciit- As the 
thick solid line in Figure 5, which is based on the value clef = —20.9 nm for the Florence experiment [27], agrees quite 
well with the crosses of the experiments, we restrict ourselves here to fit clef to the data of the Hamburg experiment. 
In order to estimate the uncertainty of the latter, we have plotted in Figure 8 a dot-dashed (dashed) line comprising 
all crosses (dots) as a upper (lower) limit. The mean value for qbf is chosen central between both limits and is 
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represented by the solid line in Figure 8, which roughly separates the dots and the crosses. This estimation method 
is justified since those stable (unstable) points, which are located above (below) the solid line, extend with their total 
error AN = ^/ANg + ANp over the solid line to the stable (unstable) sector. The new value a B F = (-16.82 ± 030) 
nm differs only 12% from the old value <jbf = —15.0 ran, whereas the critical particle numbers between the dots 
and crosses are by a factor of around 3 smaller than those of the old thick solid line in Figure 4. Note that this 
new value for the interspecies s-wave scattering length leads also to a new value of the minimal fermionic particle 
number 7V^ m i n = 1, 95 • 10 5 below which the mixture is stable in accordance with the experimental data of Figure 8. 
Furthermore, we remark that the old value of cibf was determined by the Sengstock group with the help of Eq. (55). 
However, instead of evaluating the respective fermion number Np by integrating the fermionic density (35), they used 
the ad-hoc approximation [43] 

(0) 9BF (0) 

HF = V F Mb, ( 58 ) 

9BB 

where the noninteracting chemical potentials /J^ and /JjP are related to the respective particle numbers by Eq. (39). 
The above relation between disturbed and undisturbed chemical potentials is not correct. 



V. CONCLUSION 



Applying the Thomas-Fermi approximation in order to solve the Gross-Pitaevskii equation (17) for a stationary 
Bose-Fermi mixture in Section III reduces it from a differential to an algebraic equation (34). This allows to obtain an 
analytic expression for the bosonic and fermionic density profiles. The strong attraction between bosons and fermions 
gives rise to an increase of the particle densities within their overlap, accompanied by a shrinkage of the BEC and 
the Fermi gas cloud. In order to test the validity of the Thomas-Fermi approximation, we have plotted the kinetic 
energy of a boson together with its potential energy due to the trap, its intraspecies and its interspecies interaction 
energy, and the chemical potential as the total energy of a boson. This reveals the Thomas-Fermi approximation to 
be very good over a wide bulk range of the condensate cloud. The kinetic energy plays a significant role only in the 
outermost 10% of the BEC cloud. 

Furthermore, we have found that the particle densities for both species become complex at sufficiently large particle 
numbers of bosons and fermions. We interpret this as a loss of the stability against collapse. The imaginary part of 
the density can be regarded as the decay rate of the described species and starts emerging in the trap center, where 
the densities have their maximum. 

Beside the stability diagram arising from the complex solution of the Gross-Pitaevskii equation in the Thomas- Fermi 
approximation, we have evaluated the stability border within a variational method by extremizing the grand-canonical 
free energy (29) for a 87 Rb- 40 K mixture with a Gaussian test function. The resulting lines for the variational method 
with a ratio of the oscillator lengths according to the Thomas-Fermi approximation show the same behavior as the 
ones of the Gross-Pitaevskii equation. Both lines are located close in the iVs-iVV-plane and agree well. Finally, by 
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comparing the calculated stability borders with the experimental values, we have found for the Florence experiment 
that the stability border of the variational method and of the Gross-Pitaevskii equation is in good agreement with the 
experimental results. For the Hamburg experiment, however, we have obtained a discrepancy between the calculated 
lines and the experimental results to which we have fitted the interspecies s-wave scattering length to clef = (— 16.82± 
0.30) nm. Despite this, there remains a discrepancy to the s-wave scattering length qbf = (—20.9 ± 0.8) nm of 
the Florence experiment. As both experiments deal with a 87 Rb- 40 K mixture, the interspecies s-wave scattering 
lengths should coincide. A possible explanation could be that the mean-field theory in this paper is developed for 
a stationary Bose-Fermi mixture and does not include dynamical aspects of the mixture. The dynamical behavior 
becomes important for a rapid cooling of the mixture which leads to a fast increase of the condensate density due to 
the occupation of all bosons into the ground state in a short time. We conclude with the remark that such dynamical 
effects within the Bosc- Fermi mixture could be further investigated with the help of the coupled equations of motion 
(13) and (15). 
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